Correlation functions for a two-dimensional electron system 
with bosonic interactions and a square Fermi surface 



0^ 



O 
(N 



00 

o 

On 



O 
O 



X 



3.0. Fjaerestad/ A. Sudb0^ and A. Luther^ 
Department of Physics, Norwegian University of Science and Technology, N-7034 Trondheim, Norway. 
Nordisk Institut for Teoretisk Fysik (NORDITA), Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark. 

(February 1, 2008) 



We calculate zero-temperature correlation functions for a model of 2D interacting electrons with 
short-range interactions and a square Fermi surface. The model was arrived at by mapping electronic 
states near a square Fermi surface with Hubbard-like interactions onto one-dimensional quantum 
chains, retaining terms which can be written in terms of bosonic density operators. Interactions 
^ I between orthogonal chains, corresponding to orthogonal faces of the square Fermi surface, are ne- 

■ glected. The correlation functions become sums of Luttinger-type correlation functions due to the 

bosonic model. However, the correlation function exponents differ in form from those of the Lut- 
tinger model. As a consequence, the simple scaling relations found to exist between the Luttinger 
model exponents, do not carry over to the leading exponents of our model. We find that for repulsive 
effective interactions, charge- density wave/spin-density wave instabilities are dominant. We do not 
consider d-wave instabilities here. 



I. INTRODUCTION 



(— I In recent years, research in theoretical condensed matter physics has focused intensely on the possible breakdown of 

Q Landau's Fermi-Hquid theory in strongly correlated electron systems with spatial dimension larger than one. This has 
O ' been motivated by Anderson's observationEl that the normal state of tlae high- Tc superconductors has characteristics 
which seem irreconcilable with the predictions of Fermi liquid theoryE It is well known that one-dimensional (ID) 
systeiaa of strongly correlated electrons with gapless excitations belong to the socalled Luttinger liquid universality 
^ , class,cfQ but an analogous rigorous theoretical foundation for non-Fermi liquid physics in two-dimensional (2D) systems 
\^ ' is still lacking. 

On a number of authors have treated the 2D interacting electron gas with a circular Fermi surface using bosoniza-. 

tion and/or renormalization-group techniques .ETEI Rcsummations of perturbation tbcory have also been carried out,El 
including the case of a square Fermi surface within the parquet-approximation.EI It is invariably found that the 
Luttinger-liquid fixed point is unstable in more than one spatial dimension. The underlying physics is that, due to 
Q«,^ ' the restricted kinematics on the Fermi surface in one spatial dimension, forward scattering is singular. In highjes. 
-j—; dimensions, small-angle scattering appears to avoid such a singularity, unless the interactions are very lQtie;-ranged.ll3 
On the other hand, for a nested or nearly-nested Fermi surface such as is observed in the high- Tc cupratesjiil the phase- 
^ , space for small-angle scattering is suppressed. Under sifci circumstances, it appears that a breakdown of Fermi- liquid 
I ' theory may be feasible even with regular interactions. E£l It should be noted that approaches based on, for instance, 
■ O \ resummation of perturbation series, such as the parquet-approximation of Ref. ^, are expected to be reliable for not 
^ ■ too strong correlations. Alternatively, one could consider a regime of interaction-parameters where correlation effects 
give rise to qualitatively new phenomena that should be incorporated even in a zeroth order approximation to the 
problem. This would then complement, for instance, the work of Ref. ^. 

Recently, a 2D system of electrons on a square lattice with nearest neighbour hopping and Hubbard-like interactions 
was considered, when the Fermi surface is a square even away from half-filling.EJ The necessary requirement for this 
5J] ' to be a valid starting point is that the interactions of the underlying lattice model must be strong enough. This 
\ approach therefore does not have a free-fermion limit; the zeroth order Hamiltonian to the problem is rather that of 
free bosons. In Ref. [ij, the problem was mapped onto two sets of onc-dimensional chains, one set for each axis of 
the square Fermi surface. The kinetic energy and the interaction were separated into terms which could be written 
in terms of bosonic density operators for the chains, and terms of non-bosonic form. The bosonic Hamiltonian was 
then diagonalized, assuming no interactions between orthogonal chains, and the single particle Green's functions 
were evaluated, with a resulting Luttinger-liquid form. The non-bosonic terms in the Hubbard-like interactions were 
classified as relevant or irrelevant with respect to the solution of the bosonic Hamiltonian. Discarding irrelevant 
operators, the Hamiltonian for a particular case of sufficiently strong repulsive interactions was solved, treating the 
relevant operators as perturbations on the bosonic solution. 
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The A''-component one-dimensional electron gas has very recently been treated using bosonizationjli obtaining 
the operator dimensions for a number of non-bosonic terms that could appear in the HamiltonianO One particular 
realization of such a system is the iV-chain Hubbard model. Given the fact that we neglect interactions between 
orthogonal chains in our model, the system we effectively study has much in common with that of Ref. However, 
our work differs in one important respect, namely that we consider correlation functions of fermion operators that are 
non-local in chain indices. 

The minimal microscopic model appropriate for the description of high-Tc superconductivity is ppt-vet agreed 
upon. Much work has been done for Hubbard-like models, also for the square Fermi surface case.E£rtj Recently, 
claims were however made that non-existence of d3.2_j,2-wave superconductivity in the 2D Hubbard rnadel had been 
rigorously proved at any finite temperature, and under certain circumstances even at zero temperatureHj If that were 
to be true, the single-band Hubbard model is not likely to be a fruitful microscopic starting point for explaining the 
phenomenology of the high-Tc cuprates, where there is now broad consensus on a dominant da,2_j,2-component of the 
order parameter. In this situation, the examination of the properties of various hypothetical effective theories may 
give useful information about what directions one should take in solving this exceedingly difficult problem. 

In this paper, we investigate the ground state properties of the bosonic Hamiltonian in more detail, through the 
calculation and analysis of explicit expressions for various ground state correlation functions, which enables us to 
identify the dominant divergent susceptibilities. We calculate the correlation functions for 2kF charge-density wave 
(CDW), 2kF spin-density wave (SDW), s-wave singlet (SS), and triplet (TS) superconductivity fluctuations. The 
similarities to the results obtained for the Luttinger model are: (i) the separate terms making up the correlation 
functions have Luttinger-model form, due to the bosonic form of our model, (ii) there is a degeneracy between CDW 
and SDW fluctuations, and between SS and TS fluctuations, (iii) for effective repulsive interactions, the CDW/SDW 
instabilities are dominant. However, because the expressions for our correlation function exponents differ in form from 
those of the Luttinger model, the simple scaling relations valid for the exponents of the Luttinger model do not hold 
for our model. 

This paper is organized as follows. In Sec. |l|, we describe the fundamentals of the model under consideration. In 
Sec. |h, various ground state correlation functions for the model are calculated, and their exponents are discussed. 
Sec. [V summarizes the obtained results. 



II. THE MODEL 

In this section, we give an introduction to the model under consideration, and establish the notation that will be 
used in the rest of the paper. For a more detailed discussion we refer the reader to Ref. 

A. The kinetic energy. Mapping to chains 

We consider a 2D electron system on a square lattice with lattice constant a and N lattice points in each direction 
(although we will not indicate it explicitly, the thermodynamic limit N oo will always be implied at the end of our 
calculations). The kinetic energy operator is 

where and Cis are creation and annihilation operators for an electron with spin s at site i. The sums over i and j 
are restricted to nearest neighbours, where t is the hopping matrix element. We introduce coordinates {x,y) in real 
space and (kx, ky) in reciprocal space, where the coordinate axes are rotated -1-45 degrees with respect to the primitive 
vectors connecting nearest-neighbour sites. The kinetic energy operator can then be written 

-ffkin = -4t^cos(fca;a/\/2)cos(fc2;a/\/2)ct^cg^, (2) 

ks 

where the k summation is over the 1st Brillouin zone. For a half-filled band, the Fermi surface is a square, with faces 
defined by kx or ky equal to ^kp where kp = 7r/-\/2a, giving a vanishing Fermi energy. 

Since we are interested in the low-energy physics, the kinetic energy is linearized in the direction perpendicular to 
each of the four faces of the square Fermi surface. Consider the two parallel faces defined by kx = akp, —kp < ky < kp, 
where a = ±1. These faces will from now on be referred to as the '-f ' and '-' faces, respectively. The linearized kinetic 
energy dispersion near face a is 
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(3) 



where vpiky) = 2\plta cofi{kya / \/2) . This form is reminiscent of the dispersion for the Luttinger model, except 
that in our case, the dispersion is only valid within a finite region {—kp < ky < kp), and the Fermi velocity vp{ky) 
depends on the momentum parallel to the face, vanishing on the edges of the face. We also define the field operator 
for excitations near face a as 

Mx, y^^J^Y. (4) 

where the summation over ky is between —kp and kp, and the summation over k^ should be cut off for ak^ <ti kp. 
Similar considerations and definitions can be made for the two other faces of the square Fermi surface. 

We now consider a two-point function like {4>'' {x,y,t)il! (0,0,0)) , where the average is taken with respect to the 
ground state of the kinetic energy operator. The total field operator is written as a sum of the four face field 
operators of the type given in Eq. making the two-point function above a sum of four terms, one from each face. 
The contribution from face a will be {tpl^{x,y,t)'ilja{Q, 0,0)). It can be shown that the result for this correlator is 
exactly reproduced by the following mapping of the field operator, 

I 2N-1 

tpasix,y) \ T- y2 9(1 - l')i-'aVs{x) (5) 
V l'=0 

with the resulting mapping of the unperturbed Hamiltonian, 



EE 

s 1,1' 



ti^i, J dxtpli^ix) (^j-^ - kp^ ijai's{x), (6) 



where 



N/2-1 

= ^ E Mky)e''^^'-'\ (7) 

ny = -N/2 

Here the integer I = ^/2y/a G [0, 2N — 1], ky = 2'Knyj \/2N a, Py = irny/N, so that kyy = pyl. Furthermore, 



q(1) = — V e'Py' = H)' 1 -(-!)' 
Jy'-) OAT ^ oAr 1 ^i-^i/N 



2N ^ 2N l-e^^^/^ 

ny=-N/2 

satisiying J^i^ ~ h)g{h — I') = g{l — I')- Letting N ^ oo for fixed I, we obtain 

9(1) = ^ £l If e-^(^) ^ \6,o + ^^(1 - ^.o). (9) 

The operators ^Jais{x) describe right-moving (a — 1) and left-moving (a — —1) electrons living on fictitious one- 
dimensional chains (labeled by I) of length L = Na. In terms of operators creating electrons with definite momenta, 
these chain field operators are given by 

^alsix) = 4f E CaUka;), (10) 

fcx 

where the normalization indicates that x is now taken to be a continuous variable. The discrete wavevectors are 
given by kx = 2Trnx/L, with integer Ux. Note that since g{l — I') is generally non-zero for I ^ I', the mapping (^) is 
non-local. 

We take the non-interacting system at half-filling with linearized dispersion as our unperturbed problem. Pertur- 
bations, in the form of interactions or operators taking the system away from half-filling, can in principle be treated 
within the framework of many-body perturbation theory, where their effect is expressed in terms of integrals over 
time-ordered multiple-point correlation functions. Using Wick's theorem, such correlation functions can be expressed 
as sums of products of two-point functions. Since the latter are correctly reproduced by the mapping (||) (with (||) 
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governing the dynamics), the entire perturbation series will be correctly generated. This, together with the fact that 
standard ID bosonization can be applied to the chains, is what makes the mapping useful. 

Electron creation and annihilation operators belonging to different faces and/or chains anticommute, by definition. 
Thus the associated density operators Pais{nx) = J2k' '■ '^aisi^'x + kx)cais{k'^) ■ {■ ■■■ '■ means normal-ordering) 
commute. For the '+/— ' faces this is expressed by the commutator 

[Palsinx), Pa'l's'i-n'^)] = -arixSa^a'SlJiSs^s'Sn^^n'^. (11) 

This is simply the result for the usual single-chain case generalized to the case of many chains, the only effect of the 
generalization being the introduction of the factor Si jt . The Fourier transformed operators 

2N-1 

Pas{n) = Paisinx)e~'P''\ (12) 
where n — {nx,ny), then obey the commutation relations 

[Pasin), Pa's'i~n')] = -2Nanx5aa'Sss'5i{M'- (13) 

It was shown in Ref. |lj that in the presence of interactions, the terms with I' I in (||) are irrelevant in the RG 
sense. Therefore we only keep the diagonal part, which we define as H^^^ ^. This will be a sum of independent chains. 
Bosonization is performed for each chain, giving 

Hkin,a = ■ Pepsin) Pas{-n) (14) 

where vq = tu = 2\/2ta/TT, and the sum over riy is between —N and iV — 1 (this will be the summation range for Uy 
in the rest of the article as well). 

Completely analogous considerations can be made for the two faces with ky = ±kp, the only difference being that 
the associated chains will then lie along the y direction, taken to be continuous, while the x direction is kept discrete. 
Boson density operators for perpendicular faces commute. 



B. Interactions. Solution of the bosonic Hamiltonian 



The form of the interactions which are included in the bosonic Hamiltonian can be extracted from the Hubbard 
model interaction, given by U J^i where U is the on-site repulsion and n^g is the number operator at site i with 

spin s. Writing the field operator as a sum of the field operators for the four faces of the square Fermi surface, the 
Hubbard interaction will contain terms which couple perpendicular faces and terms which only couple parallel faces. 
The former terms are not included in the bosonic Hamiltonian, while from the latter terms, which for the ' -I- /— ' 
faces can be written 

^(^) E E Eff(^-^i)3(^-^2).9(/-/3)5(^-^4) :V'l,,,t(a;)V'a.i.T(^)^l3;3i(^)^"4Ui(a;) (15) 

x,y °'l<^2 lil2 

we retain the bosonic contributions, which have the form of density-density interactions, obtained when (ai,^i) = 
(0:2,^2) and (0^3,^3) — (a4,Z4). (Note that since perpendicular faces are decoupled in the bosonic Hamiltonian, it is 
from now on sufficient to focus on one of the two sets of parallel faces, and we choose the ' + /— ' faces.) The bosonic 
contributions can be written 

^EE^^^"!/) ■ ^"T(")P"'i(-") (16) 

aa' n 

where f{ny) = (1— \ny\/N)/2 is the Fourier transform oi g^{l). Introducing charge (z/ = p) and spin (1/ = cr) operators 
as {r]p = 1, = -1) 

V2va{n) = Pa]{n) -f ri„pai{n), (17) 
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with commutation relations (i/, v' = p, a) 

[va{n),i''a'{-n')] = -2Nmix5yui5aa'Sii,n', (18) 
the most general interaction of form (^6|) can be written on bosonized form as follows: 

hb^^Y. EE/'K)^-f^-'("-) ■■ ^^i^wi-^) ■■ (19) 

u—p,a aa' n 

where we have allowed for different coupling functions uj^^,{nx) for the different scattering processes. Hence, this 
bosonic interaction is a bit more general than that extracted directly from the Hubbard model. In analogy with the 
ID case, a cutoff on the momentum transfer in the x direction will be needed in order to get a well-defined solution; 
this makes the couplings n^j-dependent. Symmetry considerations dictate that there can be at most four different 
coupling functions: u'fl, u[^l, U^^}, and U^l. 

The kinetic energy Hq = X)q -^kin a '-■^^ written in terms of charge and spin operators, 

^0 = 1^ E EE ^ ^a(")^.(-n) : . (20) 

The Hamiltonian Hq + Hb, which shows spin-charge separation, will be referred to as the bosonic Hamiltonian. It 
may be diagonalized separately in the spin and charge sectors, using the canonical transformation exp (X^i/ S^){Ho + 
Hb) exp (- S^) = H^, where 



^^=2^E^'^+(-)-(-")' (21) 



and £,v{n) is chosen so that the off-diagonal terms in the transformed Hamiltonian Hd vanish. This gives 

^^ = ^EEE"-(") ■^"(")^"(-") (22) 



2NL 

u a n 



where the parameters S,v{n) and the renormalized velocities v^{n) are given by 
exp [2^1. (n)] 



1 + ar^^U^^in,) + u'^+l{n,)]P{ny) /nv^ ^^^^ 



\ 1 + arj^[U^;l{n,) - U^+l{n^)]P{ny) /t:vq ' 
^ = \j {l + aT^,U^4{n,)Piny)/nvoy - (ary^C^t^ (".)/' (^,)At;o)'. (24) 

Terms in the microscopic Hamiltonian which are not of the type included in Ho will, if they are irrelevant, only 
lead to modifications in the numerical values of the coupling functions, leaving the structure of the low-energy theory 
unaffected. In this case, it is these 'effective' function values that should enter in the formulae above. As already 
pointed out below Ku. (p^, the non-diagonal terms (i.e. / ^ I') in the kinetic energy (^ are irrelevant perturbations 
with respect to Hd-^ We emphasize that for the-case of the 2D repulsive Hubbard model, Umklapp scattering, which 
is of non-bosonic form, is a relevant interaction.E3 The general philosophy of the approach presented here, is to view 
Hb) as a zero-order Hamiltonian for a perturbation treatment of the relevant non-bosonic interactions generated by 
the underlying microscopic lattice Hamiltonian (which may be different from the Hubbard model). 

From the high- Tc point of view, the interesting case is doping away from half-filling. In Ref. 12. it was shown that 
for sufficiently strong interactions, the only effect of the operator causing deviations from half-filling was to shift the 
value of the Fermi wave vector, not changing the square shape of the Fermi surface. In the remainder of the paper, 
we will assume that we are in this strong-interaction regime with a square Fermi surface. 

III. CORRELATION FUNCTIONS 

We now turn to the calculation of correlation functions for various fluctuations of interest in the ground state |) of 
the bosonic Hamiltonian Hq + Hb- The dominant instability is identified by the correlation function with the slowest 
space-time decay, thus having a generalized susceptibility with the strongest divergence as T ^ 0. 



5 



A. Definitions of correlation functions. Calculation of their y-dependence 



The generic form of the correlation functions we consider is 

R>{r,t)^ ~t{\0{r,t)0^ {0,0)\), (25) 

where 0(r, t) is the operator for the fluctuation under consideration. We will consider the following fluctuations: 
2kp charge-density wave (CDW), 2kF spin-density wave (SDW), s-wave singlet superconductivity (SS), and triplet 
superconductivity (TS). The SDW operator has 3 spatial components, and the TS operator has 3 components corre- 
sponing to the total spin in the z direction. Assuming spin-rotation invariance, the 3 associated correlation functions 
should be equal, although the expressions will be formally different due to the abelian bosomjation used here.Q The 
definitions of the operators O are simple generalizations of the corresponding ID deflnitionsO'cl 

OcDw{r^ = Y.^+A^^-A^, (26) 

S 
S 

OsDwA^ = E^^L(^>-^-(^^)' (28) 

S 

Ossirl = ^ 5]^^,,(f)V'L,_,(r-), (29) 

^ S 

OTS.+i{r) = i^lA^)^lA'^- (30) 

We have included the definitions of two of the components of the SDW operator, since the equality of the associated 
correlation functions will give a condition for determining the parameter S,a{n). 

Introducing the mapping (|^), and observing that the bosonic Hamiltonian conserves the number of electrons with 
a given spin on a given face and a given chain, it is clear that each non-zero term in the correlation functions has the 
form (discarding prefactors) 

Fi3,s,s, {h -l2,x,t) = (IVl,,,,,, {x, i)V^^,,,,,, (x, t)iri^s, (0, 0)^+,h,s^ (0, 0)1), (31) 

where /3 = -1-1 (—1) indicates an annihilation (creation) operator. Here we have anticipated that -F is a function only 
of the product of the spins, and of the separation between the chains. In terms of F, the correlation functions are 
given by 

2N-1 2N-1 

i?>^^(r,i)cx ^ Hi{hM)F++{h-l2,x,t), (32) 

il=0 (2=0 

2N-1 2N-1 

RsDwJ^^t)^Yl E ^'('l''2)^+-('l-'2,X,t), (33) 
2N-1 2N-1 

RsDwJr,t)^Yl E ^'('i''2)^++(^i-^2,^,i), (34) 

il=0 (2=0 

2N-1 2N-1 

Rssi^,t)^Yl I] iri('i,'2)F^— (?i-?2,a;,i), (35) 

il=0 l2=0 

2N-1 2N-1 

^k+i(^'*)« E E ^'('i''2)^-+('i-'2,a:,t), (36) 

il=0 l2=0 

where Hi{li,l2) = g{—li)g{l ~ h)g{—h)9{l — h)- It is seen already at this point that in this model, the CDW and 
SDW correlation functions are equal. 

In (^2|)-(^^ we first sum over all terms with a fixed value of /q = /i — 12, and then sum over the appropriate values of 
Iq. Using Hi{li,l2) — Hi{l2,li) and F(/o, x, t) — F{—Iq, x, t) (the latter property is established in the next subsection) 
we obtain 
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2N-1 



R{r,t) CKK{l,0)F{0,x,t) + 2 ^ K{l,lo)F{lo, x,t), 



ln = l 



where we have defined 



2N-1 



K{lJo) = J2 Hiih,h-lo) 



(37) 



(38) 



h—to 

For /o > 1 it is convenient to write K{1, Iq) — K^{1, l^) — K_{1, Iq), where 

2N-1 



ii=0 

i^_(?,?o)-^if(('i,^i-^o). 



(39) 



(40) 



In the calculation of i4r+, the summation over li must be done before the N oo limit is taken. This is because the 
expressions (^) are derived by assuming that I is fixed, and therefore finite, so that l/N — s- as A'^ — s- cx). They are 
therefore not correct when I is of order TV, and such values of li are indeed included in the sum in K^. 
The calculations of and K- are rather lengthy, so we only give the final results {5i,ii = 1 — di^i): 



1 - 1 _ 1 _ (-1)'+^ 



(41) 



K-{l,lo) 



sin(7r/o/2) 



27rZo 



cos2(7r?/2)cos2(7r/o/2) 



"2(V'( 



V,(i±i)-^(i±^) + V.(i) 1 fibi^) , V'(^)~ 



lo{l + lo){l - lo) 



2//q \ I -\- Iq 



l-lo 



5lo 



l+lo 
2 



/3 



4/2 



+ S, 



Llo 



^ + ^(i)_^(i^)„;V,^(i^) 
4/3 



(42) 



Here ipiz) is the digamma function, defined as '0(z) = T'{z)/T{z), where T{z) is the gamma function.Ell For large 
arguments, ip{z) ln(z), which in this context must be regarded as a constant, since it does not contribute a power- 
law behaviour with non-zero exponent. The asymptotic behaviour of K{1, Iq) is therefore to leading order given by 
K{l,lo) oc 



B. The X and t-dependence of the correlation functions 

The explicit calculation of F is performed by bosonization of the chain field operators ,01 writing 

i;^is{x, t) = 1 [/t^^e-"^°-(^'*)+^"'=-^ (43) 
V27re 

where e is a short-distance cutoff which should be sent to zero at the end of all calculations. In principle, the phase 
field Lp contains both finite modes (rij. ^ 0) and zero modes (n^, = 0, corresponding to number operators), but the 
zero modes may be neglected in the thermodynamic limit.a Furthermore, it is easily seen that the ladder operators 
U do not give rise to any minus signs in the expectation values we consider, because there is no need to reorder the 
C7's before eliminating them using UaisUl^is — Ul^igUais — 1, and they can therefore be neglected. The phase field is 
ipaisix, t) = ^aisi^^ where we have defined 

^:,,{x,t) = ^^£:^.„(n,Oe-'=--^«\ (44) 
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where Hps = 1, k^-s = s with s = 1 (—1) for spin up (down). The time dependence is given by Va{n,t) = 
^i{Ho+HB)t^^(jl-^^'i(Ho+HB)t _ Ymia i^l{n) Va{~-n), we obtain ip^^i^{x,t) = -ipais{x,t) . This gives 

Fp,s,sAh-h,x,t)^{2TTe)-^exp[-ikFx{l+f3)] Y[ exp[E^p^,^,^{h - l2,x,t)], (45) 

l^ — p.(T 



Next, we perform the canonical transformation with the operator (|21|). Its effect on the spin and charge operators is 
given by 

e^-i/„(n,i)e-^- = ^ K\n)Mn)e'^'^^^"(''^\ (47) 

A=±l 

where 

/i^'"(n) = coshf^(n), K-^'^in) = -sinhi^{n). (48) 
The canonical transformation will mix the operators for the '+' and '-' faces in the phase fields: 

e«''^:;,,(x,i)e-^^= E ^lU^.t), (49) 



where 



A=±l 



^"^isM = -^^f:^/,^^(n)z.,(n)e^^'=-'^(")*-'=--^"''. (50) 

n 

This gives 

eMEl,^,^{h-h^x^t)]^ n (G|e^+'i=i(^'*)e''^-'2=2(-*)e-''^-2e2(°.«)e-*^^^ (51) 

A=±l 

where \G) = exp Si,)\) is the ground state of Hd- 

The calculation of the expectation value is done by writing the ^p^ig fields as a sum of two contributions, one over 
positive Ux and one over negative Ux- One of the contributions will contain creation operators, the other will contain 
annihilation operators. By using the relations 

gAgB ^ ^A+B^[A,B]/2^ (53) 

(53) 

we can move the exponentials containing creation operators to the left and the exponentials containing annihilation 
operators to the right. In doing so, we will generate a string of c-numbers due to the reordering of the exponentials 
which we can take outside the expectation value brackets. Since the annihilation operators acting to the right and 
the creation operators acting to the left annihilate the ground state, the expectation value of the operators which are 
left inside the expectation value brackets is just 1. 

First, the phase fields are split into the two different contributions, 

V^isM = ^S;^(x,t) + (54) 
where the superscripts + and — denote creation and annihilation, respectively, and 

(p'2if{x,t) = ±^^^'£l!^/i^^(n)zyA(±An)e±[^'=^"-(")*-*^'=-^-^^P'''l, (55) 
2v2A' ^ T^x 

n 

where we also used the evenness of S,v{n), and where the prime after the summation sign indicates that only positive 
Ux should be summed over. Then we use ( p2|) to write each exponential in ( pT|) as a product of exponentials of creation 
and annihilation parts of the phase field. For this, and further calculations, we will need the commutator 
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fl y^'f l/jOA/^N^a'A/^x ifc,i;„(n)(t-t')-iA[fc,(x-2:')+p„(i-i')]_ 



4iV ^ 

n 

Using ( ^2|) on all exponentials in (|T]) we generate the factor exp [C] where 

1 



■ cosh 2^y(n), 



(56) 



(57) 



where we also used — 1. Using ( ^3| ) to move all annihilation operators to the right in (|5l|), we generate the factor 
exp[D^;\^^^(?i -h,x,t)\, where 



4A^ ^ 

n 

sinh2C,(n)(l - e-''^'=='^''(")*+*^'=-^) cosp.y(Zi - ^2) + cosh2C.(n)e'^'=-''-(")*+'^'===n 

The total exponent is 

1 

A=±l 



(58) 



1 



[cosh2^y(n) - PKysiKys2 sinh2^i,(n) cospj^(/i - /2)]- 



kx 



(> 



-ik^Viy {n)t 



COS k'^oc 



(59) 



A momentum transfer cutoff must be imposed on the interactions. For fixed Uy, the treatment is as for the ID case. 
In the equation above, we express the integral as a sum of a non-interacting and an interacting part, by writing 
[•••] = 1 + ([•••] — 1). In both parts, we then multiply the integrand by [e"'^-' + (1 — e~^''("f)'^==)] to obtain two 
terms for each of the two parts. (For simplicity, we take the cutoff function A^(ny) to be the same for both parts; since 
the non-interacting part is independent of /?, S1S2 and li — I2, this implies that the cutoff function be independent 
of these variables.) In the first term, the main contributions will come from kx ^ 1/A^(ny), where we can replace 
iyvin) — > S,u(nx = 0,ny) = S,u{ny) and Vy{n) v^irix = 0,ny) = Vy{ny). In the second term, which can be neglected 
for the interacting part, the main contributions come from ^ 1/A^(nj^), so we can replace Vi,{'n) — > vq there. This 
gives 



Ef)^,,,^ih - l2,X,t) 



4- 



(e + ivoty + x^ 



^ I dpy ^ In 



{Ai,{ny) + ivi,{ny)ty + x^ 



(A^(n^) + iVy{ny)tY + -t^ 



where we have defined 

^/3,siS2("y' ^1 - ^2) = -1 + cosh2^i.(nj^) - PKysiKus2 sinh2fi,(nj^) cospyih - h)- 



(60) 



(61) 



Introducing nj,-averaged velocities Vi, and cutoffs A^, and adding and subtracting terms in the integrand, we may 
write 



n exp[E^p^^^,^ih-h,x,t)] 



l/—p,(T 



n 



{e + ivot)^+x^ \{Ai, + ivyt)^ + X 



v—p.a 



(A, + tvot)^ + x^ ^ 



A2 



exp 



/ dpy <iln 



(A^ + iv,A)^ + x^ 
{A^{ny) -^iv^{ny)tY +2;^] / [ (A^, + iu^t)^ + 



+ ^^.sisaK'^l - ^2) In 



{Ay(ny) + ivy{ny)tf + x"^ 



where we have defined the charge and spin exponents 



_{Ay + ivot)^ + x^ 
(A^ + w^t)2 



A2 



1 r 

^~ Jo 



(62) 



(63) 



Summing up the results so far, the correlation functions (p2h-(M) are given by Eqs. ([37])- (42), (}45|), and (|6l|)-(|6 
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C. Correlation function exponents 



From Eqs. ( P3|) and ( |34| ) we see that the correlation functions for the x- and z-components of the SDW operator differ 
only in the sign of siS2- Assuming spin-rotation invariance of the underlying model, these correlation functions must 
be identical. From Eq. (|6ll ) it is seen that this is obtained if £,a{ny) = 0, which gives ^^^^(ny, Iq) = d'p ,._^g^{lo) = 0. 
Consequently, the correlation functions become independent of S1S2, and the S1S2 subscript may therefore be omitted 
from now on. The s pinrr ptation invariance also has the effect of making the SS and TS correlation functions identical. 
As for the ID case,cSE3 a more sophisticated (renormalization-group) treatment is needed to lift the degeneracy 
between CDW and SDW fluctuations, and between SS and TS fluctuations. 

For i = it is seen that the factor exp[ • • • ] on the rhs in ( |6^ ) is unity for v — a, and for = p its leading 
behaviour for x — > 00 is independent of x. Thus the exponent of the leading a;-dependence of ( |6^ ) is — 2[1 + d^(/o)]. 
For asymptotic values of t the analysis is not so straightforward, since both the square roots and the factors exp[ • • • ] 
will depend on x and t to leading order. However, as long as all velocities involved are non-zero (which we assume to 
be the case) , it seems reasonable to neglect this additional dependence, so that the following leading-order asymptotic 
approximation may be used: 

''''(^°'"'')"(7TWT^l,(A;T7#T^j ■ 

We see that Fp{lo,x,t) has a form that resembles that of the correlation functions of the ID Luttinger model. 
However, the exponents d^(^o) depend on the chain difference Iq, and in general have a rather different form than 
the Luttinger model exponents. The exceptions are the exponents for the equal-chain terms Iq = 0. We have, in our 



CDW/SDW: d^^^^ = -i(l-e-2«.) (ID), ^^^^^(^0 = 0) = /; dp, (1 - e-^C^K)) (2D) 
SS/TS: rfP_^ = -i(l-e+2«.) (ID), dP_^(/o = 0) = -^L /; rf^^ (1 _ e+^U".)) (2D) 



(65) 



It is seen that the form of these exponents is very similar in the ID and 2D case, the only difference being the averaging 
over py in the 2D expressions. If this had been the set of exponents which determined the leading behaviour of our 
correlation functions, then the analogy to the Luttinger model would have been very strong indeed. However, things 
are not quite that simple. Each correlation function is given as a sum over contributions from different Iq- Since for 
all Iq, K{1,Iq) oc 1/P for asymptotically large I, the asymptotic behaviour of the various contributions differ only in 
their a;-i-dependence. Thus, for a given correlation function (i.e. a given value of /3) the leading asymptotic behaviour 
comes from the value of Iq which gives the smallest fi^(Zo), i-e. we must minimize 

p r 

~ TT I '^Py smh2^p{ny)cospylo (66) 
Jo 

with respect to the integer ^q. Assuming repulsive effective interactions (i.e. £,p(ny) > 0), it is indeed true that ^^(Zo) 
is minimal for the term with Iq = 0, which therefore gives the slowest decay in the CDW/SDW correlation functions. 
However, for d'Liilo) the term Iq = actually gives the largest exponent; the minimum exponent for the SS/TS 
correlation functions comes from some non-zero Iq, and is therefore not of the form in (|65|). Consequently, the relation 
between the leading exponents for the CDW/SDW and SS/TS correlation functions is not as for the ID Luttinger 
model, where these exponents are related through simple scaling relations. Since (i^(0) < d'Li{lQ), regardless of Iq, we 
conclude that the CDW/SDW fluctuations are dominant. This is the same conclusion as for the Luttinger model for 
repulsive interactions, although the detailed nature of the leading exponents is different. 

Throughout this paper we have assumed that we are in the region of parameter space where the square shape of 
the Fermi surfajsa is stable with respect to doping away from half-filling. This requires, when spin-rotation invariance 



is invoked, thaul^ 



1 



sinh ^p{ny){l — cosp^Zo) > 2 for all Iq ^ (67) 



(since = 1 then minimizes the expression on the Ihs of (|67|), an equivalent requirement is that this expression be 
> 2 for ^0 = !)• Let us define U ~ u\^l{0)a/ Attvq and 7 = u\^}_{0) /u\^l{0) . Having repulsive interactions requires 
> 0, 7 > 0. Furthermore, Eqs. ( |2^ ) and ( [24| ) are well-defined and real only when 7 < 1 -I- 1/U. In order to have 
satisfied, there is a lower bound on 7, which is easily found by assuming 7 < 1, and noting that for fixed 7, 



sm. 



1 ^piriy) is maximal when U ^ 00, giving 
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The condition ( |67| ) then requires 7 > 7inin = ~ 0.9798. There is also a 7-dependent lower bound on U . The 

lowest such bound is found numerically to occur at the parameter space boundary, for U ~ 376, corresponding to 
7 = 1 + 1/U ~ 1.00266. As 7 is decreased from this value towards 7inin, the lower bound on U is found numerically to 
increase towards infinity. A lower bound U ~ 376 corresponds to U^^{0)a/TTVQ ~ 1500, an extremely high value. We 
do not know how large the bare couplings must be in order to renormalize to such high effective values. It may be that 
inclusion of scattering between orthogonal faces could reduce the effective values needed, but that is not clear. In the 
region of parameter space where ( |67| ) is satisfied, we find that Iq = 2 minimizes d^j^{lo) if 7 < 1 (the expression ( |6^ ) 
is then negative), and Iq — 00 minimizes d^i(lo) if 7 > 1 ((|6^) is then 0). We do not have any physical explanation 
for why these particular values of Iq minimize the SS/TS exponent in each particular case. 

IV. SUMMARY AND CONCLUDING DISCUSSION 

We have evaluated various ground state correlation functions for a 2D bosonic Hamiltonian with spin-charge sepa- 
ration, introduced in Ref. |l2|. The model was arrived at by mapping the problem of 2D electrons on a square lattice 
with nearest-neighbor hopping and Hubbard-like interactions on a square Fermi surface onto two orthogonal sets of 
ID chains. In the presence of interactions, it was shown in Ref. ^that single-particle hopping between parallel chains 
was irrelevant; furthermore, for sufficiently strong interactions, the square shape of the Fermi surface was preserved 
even away from half-filling. The bosonic Hamiltonian resulted from neglecting interactions between orthogonal chains 
and parallel-chain interactions of non-bosonic form. 

In the mapping used, the physical 2D field operator is written as a sum of field operators on the ID chains. It 
is the 2D field operators which enter in the definitions of the correlation functions. As a consequence, the resulting 
correlation function is a sum of many terms, each of which takes the form of a Luttinger-model correlation function 
in the time direction and the spatial direction parallel to the chains. The chain indices enter in the prefactors and in 
the exponents of the Luttinger-model correlation functions. In order to find the leading behaviour of the correlation 
function for large times and distances, one must identify the term with the smallest decay exponent by minimizing the 
exponent with respect to the chain separation. The leading exponents were found to differ in form from the Luttinger 
model exponents due to the 2D nature of the problem. Specifically, the simple scaling relations valid for the Luttinger 
model exponents were not valid for our model. 

The fluctuations we considered were 2kp charge-density wave and spin-density wave, s-wave singlet superconduc- 
tivity, and triplet superconductivity. The condition of spin-rotation invariance fixed the spin part of the correlation 
functions, making only the charge part non-trivial. The CDW and SDW correlation functions turned out to be de- 
generate, and so did the SS and TS correlation functions. As for the Luttinger model, the CDW/SDW fluctuations 
were dominant for effective repulsive interactions. 

One may ask whether our conclusion concerning the absence of divergent superconducting fluctuations for repulsive 
interactions holds for pairing in any symmetry channel. In view of the Kohn-Luttinger effectJ23 one would expect to 
see T = superconductivity in some, albeit possibly very high, angular momentum channel. It may be that possible 
superconductivity is either hiding in higher spin-triplet channels, or in higher singlet-channels, not considered in this 
paper. Since Oss(r) is a local object, it does not contain a (ia;2_j,2 -component, and consequently our results do not 
rule out the possibility of a superconductive instability in this channel. 
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